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. Abstract. We use the Ultra-relativistic Quantum Molecular Dynamic model (UrQMD v2.2) to 

' study forward-backward fluctuations and compare our results with the data published by the PHO- 

C3 I BOS experiment. The extracted effective cluster multiplicities show a clear centrality dependence 

within the present hadron-string transport approach. This behavior is not reproduced with models 
not taking into account final state rescattering. 
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^ INTRODUCTION 

00 , 

O . One of the main goal of the relativistic heavy ion program is to understand the nature of 

hadron production mechanisms (e.g. parton coalescence, string fragmentation or cluster 
decay). Numerous data suggest the formation of a quark gluon plasma (QGP) during 
the collision of two heavy nuclei at ultra-relativistic energies. Using correlations and 
fluctuations to probe the nature and properties of the highly heated, high density matter 
^ ' created in the course of these collisions has been proposed by many authors, see for 

Ch ■ example 1 1, 2, 3., ^ ^ ^ i ^ ^ JOJ. In particular, one can get insight about cluster 

decay with the help of such fluctuation studies. 
kJ [ The UA5 experiment performed a study of cluster size in p — p collisions by analyzing 

'j^ ■ forward-backward charged particle multiplicity fluctuations lUlll . They found a cluster 

multiplicity around two, in line with the expected result in case clusters correspond to 
decaying resonances (e.g ^ + K~). A similar analysis was recently performed by 
the PHOBOS experiment for Am +Au reactions at ^/s^ = 200 GeV lfI3l . 

In Refs. [2, 3], a simple model was introduced to extract the effective cluster multi- 
plicity ^eff from the PHOBOS data. Within this approach, K^ff is found to be 2.7 for mid- 
peripheral events and ^eff ~ 2.2 for central events. The value of ^eff in central collisions 
is close to the p — p value reported by UA5 ifllll . Note that all measured cluster multi- 
plicities ^eff are larger than the one computed for a hadron resonance gas (^hg = 1-5) 
ll3] . indicating that the measured charge correlations can not be described by simple 
statistical models based on hadronic degrees of freedom. The STAR collaboration has 
measured the long range correlations and it was shown that also within a string fusions 
approach, the data for central A +A reactions can not be reproduced fl?]. 

In this study we calculate a baseline estimate for forward-backward fluctuations 
based on the microscopic hadronic transport model UrQMD. For a complete review 



of the model see fls*]. The centrality, rapidity 77 and rapidity window Arj dependence 
of the dynamical fluctuations are studied and interpreted in term of effective cluster 
multiplicities. For this analysis, 5 x 10^ pp and minimum bias Au+Au events at a/saw — 
200 GeV were used. 



FORWARD-BACKWARD FLUCTUATIONS 

In this section, we introduce the variable C that measures the asymmetry between 
the forward and backward charges. We define two symmetric rapidity regions at ±77 
with equal width Arj . The number of charged particles in the forward rapidity interval 
rj ±Ari/2 is Np while the corresponding number in the backward hemisphere —77 ± 
At] /2 is given by Nb- We define the asymmetry variable C = {Nf — Nb) / ^Nf +Nb, in 
each event. The variance of the charged particle multiplicity in the forward hemisphere 
is given by Opp = (Np) — (Np)^ and similarly for the backward hemisphere Dbb = 
(A^l) — {Nb)^. We also introduce the covariance of charged particles in both hemispheres 
by DpB = {NpNb) — {Np){Nb), where (....) stands for the average over all events. The 
PHOBOS measure of the dynamical fluctuations can be written as 

ai = (C2) - (C)2 ^ Dpp + Dbb-2Dpb 
c ^ ' \ ' (Np+Nb) 

Recently, STAR ifTill reported preliminary results of the so called correlation strength 
parameter b = DpB/Dpp. The effective cluster multiplicity ^eff is proportional to o^, 
such that if b=0, then the covariance DpB vanishes. In this case we have = K^ff. We 
emphasize that ^gff should be understood as a product of the true cluster multiplicity 
times a leakage factor ^ that takes into account the limited observation window At] 
. The event by event fluctuations of the asymmetry parameter (variance) in the 
absence of any correlations among the produced particles will be C7^ = 1 . If there is only 
long range correlation, then < C7^ < 1 . 



COMPARISON TO DATA 

We compare the UrQMD model calculations with the existing experimental data re- 
ported by PHOBOS for Au+Au collisions at = 200 GeV in Figs. [H |2l Also 
shown here are the HIJING resuk as reported by PHOBOS [12]. 

Fig. d depicts as a function of the pseudo-rapidity T] while the window is fixed at 
At] = 0.5. In Fig.|2l we show as a function of At] while T] is fixed at T] = 2. Central 
data correspond to the 0% — 20% most central events while peripheral data correspond to 
40% — 60% events as determined from the number of charged particles at mid-rapidity. 
We find that the present transport approach is able to reproduce mid-peripheral PHOBOS 
data roughly for both pseudo-rapidity T] and window At] dependence (Fig.[Tl right and 
Fig. 121 right). As a function of the pseudo-rapidity window At], UrQMD can mimic the 
mid-peripheral experimental data, however, it fails to reproduce the central one. 
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FIGURE 1. Dynamic fluctuations as a function of pseudo rapidity rj at -y/smv = 200 GeV, with Arj = 
0.5. (Left) 0% — 20% central Au+Au, (right) 40%-60 peripheral Au+Au collisions. FuU circles are 
PHOBOS 1 12] data, open squares HIJING and open circles are UrQMD results. 



As a function of rj (Fig.[T]), increases from C7^ ^ 1 to ~ 1.75 for mid-peripheral 
events. For central events, it increases from 1 to ^ 1.6. The value ^ I 
when rj =0.25 can be explained by the competition between long range and short 
range correlations which almost cancel out when the forward and backward acceptances 
are very close. The negative long range component then decreases with rj and let 
increase. 

Fig. El depicts the dependence of on the size of the rapidity window Arj. The 
experimental value of increases up to 2.8 for mid-peripheral events and reaches 
= 2.2 for central events. UrQMD result overshoots PHOBOS data for central events. 
In Ref. 10, 111], we argue that by increasing the observation window Arj one can see the 
whole cluster structure. The only process able to destroy the cluster structure in UrQMD 
is hadronic rescattering. The failure to reproduce central data can thus be seen as an 
indication for a lack of rescattering in UrQMD. 

HIJING does not yield any centrality dependence and does not correctly reproduces 
the data for peripheral events, however, surprisingly reproduces the central data. 



CENTRALITY DEPENDENCE 

In this section, we study the centrality dependence of the forward-backward fluctuations 
calculated from UrQMD for a set of Au + Au and p + p events at the highest RHIC 
energy available, as a function of rj is shown in Fig.|3](left) for different centralities. 
We observe a clear increase of from p + p {o^ = 1 . 1 at T] = 2.75) up to the 20 — 40% 
Au+Au most central events (o^ = 1.8 at 7] = 2.75). With even higher centrality, the 
behavior then changes and gets now smaller (c7^ = 1.6 at t] = 2.75 for the — 20% 
most central events). The same trend is observed as a function of the pseudo-rapidity 
window Arj shown in Fig. |3l (right), increases from p+p (a^ = 1.25 at Arj = 2) to the 
20 - 40% Au +Au most central events (a^ = 2.75 at At] = 2). In lHH we predicted that 
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FIGURE 2. Dynamic fluctuations as a function of pseudo rapidity window Arj at -^/saw = 200 GeV, 
with rj =2 (Left) 0% — 20% centi-al Au+Au, (right) 40%-60 peripheral Au+Au coUisions. Full circles are 
PHOBOS data O], open squares HIJING and open circles are UrQMD results. 




for PHOBOS data, may be reduce to 1.9 with tighter centrality cuts. This reduction 
in ~ ^eff may be regarded as an indication for cluster melting at RHIC. 

A summary of the centrality dependence is presented in Fig. |4l where is shown 
for fixed T} =2 and At] = 2. With this acceptance window, increases from = 2.1 
for peripheral events up to C7^ = 3 for the 20 — 30% most central events. Followed by 
a decrease down to = 2.7 for the — 10% most central events. From Fig. El we see 
that a model without final state rescattering like HIJING seem not to be able to mimic 
the decreasing cluster size with the most central events. 
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FIGURE 4. Dynamical fluctuations computed with UrQMD as a function of centrality , with J] = 2, 
and At] = 2. 



SUMMARY AND CONCLUSION 

In this paper we computed forward-backward fluctuations and compared UrQMD cal- 
culations results to the available experimental data measured by the PHOBOS collab- 
oration. The study of proton-proton collisions indicated that the long range correlation 
persists over a wide rapidity gap between the two rapidity hemispheres. The variance 
of the asymmetry parameter C was found to increase with increasing At] such that 
changes from o^{r\ = 2, At] = 0.25) ^ 1 to o^{r\ = 2,At] = 3) ^ 1.6, this can be due 
the saturation of the leakage factor ^ \ . 

For Au-i-Au collisions, we found that for both centrality bins 0% — 20% and 40% — 
60%, 1 for small t] . This can be seen as a cancellation between short and long 

range correlations. By increasing Tj we observed that also increases and approaches 
1.6 and 1.8 for 0% — 20% and 40% — 60%, respectively. This increase can be attributed 
to the decrease in the long range correlations. This will be true if the particle production 
mechanism does not change with T]. To see the whole cluster structure, we fixed the 
center of the observation window at 2 and allowed At] to increase. We found that 
UrQMD can reproduce the peripheral data while it overestimates the experimental 
results for central collisions. This might indicate an additional cluster melting process 
not accounted for in the present hadron-string model. 
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